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Abstract 

A perturbation scheme is discussed for the computation of the normahzation 
constant of the large order behavior arising from an ultraviolet renormalon. In 
this scheme the normalization constant is expressed in a convergent series that 
can be calculated using the ordinary perturbative expansion and the nature of the 
renormalon singularity. 
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The large order behavior in field theories arising from a renormalon is generally given 
in the form 

a„ = /s:n!n'^6o"[l + 0(l/n)] for n ^ oo. (1) 

While the constants u and Bq are calculable, the normalization constant K cannot be 
determined exactly. An infinite number of renormalon diagrams contribute to it, but it 
is not known how to sum such diagrams to all order ^, ^. 

Though the normalization constant cannot be determined exactly, we shall see that it 
can be calculated perturbatively to an arbitrary precision. The large fiavor (Nj) expansion 
is often invoked for an approximate evaluation of the normalization, and in the literature it 
is further asserted that the normalization cannot be computed without resorting to it 0]. 
However, the Nf expansion may not be considered a systematic perturbation scheme, as 
it is not proven that it gives a convergent series or even compatible with nonabelian gauge 
theories. At large Nf, the asymptotic freedom is lost, and so it is probably incompatible 
with asymptotic free gauge theories. 

We note, however, that there is a systematic method that is both compatible with 
nonabelian gauge theories and gives a convergent series. In we have discussed a scheme 
which express the normalization constant of an infrared (IR) renormalon in a convergent 
series which depends only on the strength of the renormalon singularity and the ordinary 
perturbative coefficients of the amplitude in consideration. A sample calculation in QCD 
using the radiative calculations up to three loops shows that our method gives a rather 
quickly convergent series. The purpose of this letter is to extend the scheme to the case 
of an UV renormalon. 

Let us first review the perturbation scheme briefiy in the case of an IR renormalon. 
To be specific, we consider a Green's function D{a), such as the Adler function, in QCD 
and its expansion 

oo 

D(a) = $:a„a"+^ (2) 

n=0 

The Borel transform D is then defined as follows. We first define it in the neighborhood 
of the origin as 

m = t ^b-^ (3) 

n=0 "■■ 
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and then by analytically continuing to the whole 6-plane. D{b) is known to have IR 
renormalon singularities at 6 = —n/ (3o, n = 2, 3, 4, • • • and UV renormalon singularities at 
b = n/ Pq, n = 1, 2, 3, ■ ■ ■, where Po is the first coefficient of the /3-function. Throughout 
this letter we assume there are no other singularities associated with D{b) except for those 
caused by the instantons which are irrelevant for our discussion. Now consider the first 
IR renormalon at 6 = —2/j3o. The nature of the first IR renormalon singularity is given 
in the form 

^W-^^^ for b^^, (4) 
1 + M ] Po 



2 

where u = — 2/3i//5q, and Pi is the second coefficient of the /3-function. Then the large 
order behavior caused by the renormalon is given in the form (|I|) with K = D/u\ and 
6o = —2/Pq. Note that the normalization constant becomes the residue of the singularity 
in the 6-plane. Thus we can equivalently work with the renormalon residue in order to 
study the normalization constant. 

The perturbative calculation of the residue is based on the observation that the residue 
of the closest singularity to the origin in the complex plane can be expressed in a conver- 
gent series involving only the Taylor expansion at the origin. Assuming that the strength 
of the singularity is known all one needs to calculate the residue is to move the singularity 
in consideration by conformal mapping closer to the origin than any other singularities in 
the Borel plane. 

This scheme works as follows in the case of the first IR renormalon. The closest 
singularity to the origin in the 6-plane is the UV renormalon at b = I/Pq. Using the 
conformal mapping 

Pob 
1 - Pob 

we can make the first IR renormalon the closest singularity to the origin. In the 2;-plane 
the first IR renormalon assumes 



1+1/ ^ 

D{b{z))^ for z^\. (6) 

- - z] ^ 



Now consider a function defined by 



R{z) = Dm) [ --z] . (7) 



At the IR renormalon singularity, R{z) could still be singular, but it is bounded. The 
residue is then given by 

b = R{z) \. . (8) 



.9. 

Since R{z) is analytic on the disk \z\ < 2/3 we can write in a series form by expanding 
it at 2; = 

/Ox 1+1/ 00 

(5) w 

It is easy to see that r„ depends only on with i < n, and thus calculable. 

One might question the convergence of the series (|^), since it is evaluated at the 
renormalon singularity at which R{z) could be singular. However, it should be noted 
that the finiteness of R{z) at the singularity guarantees the convergence. A numerical 
evaluation of the series using the 3-loop calculation of the Adler function in QCD shows 
a quick convergence for small Nj case For example, the first three elements of the 
series are 0.904, -0.358, 0.003 for Nf = 2, and 0.946, -0.354, -0.098 for Nj = 3. 

This example demonstrates that a renormalon residue can be expressed in a convergent 
series once the nature of the singularity is known. We now show that a similar conclusion 
can be made with a UV renormalon. In the following we shall assume that (3o < 0, for 
definiteness, and focus exclusively on the first UV renormalon since it gives the dominant 
large order behavior. The structure of the UV renormalon is a little more complicated. 
According to Parisi, it is determined by an insertion of dim=6 operators 0. To be specific, 
let us consider a Green's function A{a) and its Borel transform A{b). Generally the Borel 
transform has a branch cut beginning at the first UV renormalon on the negative real 
axis, and a quantity defined by 

ImA(a) = - lim — / dbe~^ A{b + ie) - A{b - ie) (10) 

e^O 2i J -00 L 

is nonvanishing. The Parisi ansatz states that the dominant contribution to ImA(a), for 
a — s> 0_, arises from an insertion of dim=6 operators, and is given in the form, 

M 2 

lmA{a) = Y.fi{a)Oi{a) + 0{e-^) (11) 

i=l 
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where the index i runs over all dim=6 operators Oj. And /j satisfy the renormalization 
group equation 



(12) 



where jij denotes the anomalous dimension of the associated dim=6 operators. An explicit 
implementation of the Parisi ansatz may be found in |^. Solving the RG equation one 
can write formally 



lmA{a) = - 2^ TT^^ °" IT 



1 + E 



r(z4)a 



a 
Jo 



+ 0(e 1^) (13) 



j=l - J) 

where Ki are undetermined constants while z/j depend on both f3o, (3i and the one loop 
anomalous dimension, and Cij are calculable constants depending on the higher order 
corrections on /5(a), 7(a) and Oi{a). Note that the summation within the bracket is not 
well defined; but, this point will be irrelevant in the following discussion. 
The corresponding Borel transform to ([T3|) is then given by 



M 



(14) 



in the neighborhood of the singularity at 6 = l//3o, and the corresponding large order 
behavior is given in the form 



(15) 



Our aim is to express the prefactors in a calculable, convergent series. Without 
losing generality, we may assume z/j > uj for i < j. And also, for the moment we shall 
assume all z/j > 0, and later will make a comment on the case with a negative z/j. Then 
it is straightforward to write Ki, the prefactor of the most dominant term in the large 
order behavior, in a convergent series. Since 



K, = Aib)il-Pob) 



(16) 



we can obtain a series expression for Ki by expanding the function on the r.h.s. at 6 = 0. 
The resulting series is then convergent because the function is finite at the singularity and 
there is no other singularity within the disk \b\ < —1/(3q. 



For the prefactors other than Ki we get a hnear relation among them. Using ([T^) it 
is easy to write Ki as 



Ki = [hi{b) +J2mij{b)Kj] 1^^^ 



/3o 



(17) 



where 



and 



h,ib)=Aib)il-Pobr' 



rrii 



-(1 - PobY'-'^^ [l + EtLi Cjkil - /3o&)'J , for z > j 
0, for i < j , 
with [z/j — z/j] being an integer satisfying 



(19) 



< z/j — z/j — [z/j — z/j] < 1. 



(20) 



To solve eq. (|T7|) we introduce rj(6) defined by 

r^{b) = - m,^{b)]K^ - h,{b). (21) 
j 

Note that by definition rj(6) vanishes at the renormalon singularity. From (|2T|) we obtain 

K. = E[l--^^V(/^r+rr) (22) 



where Z^^-* for a function /(6) denotes the A^-th order Taylor expansion evaluated at the 
singularity, i.e, 

(23) 



•' ^ rlhn 



n=0 



c?6" n! 



Note that /ij^'' is calculable in terms of the perturbative coefficients of A{b). From the 
definition (^) we have 



m. 



(N) 



and from (|T4D, (|TD 



r 



(TV) 



0, for i < j 

i—l oo 

E E 



'^{uj-Vi-k+l) 



for 2 > J 



(24) 



\ , r 1 , -, r(z/,- — Z/j — + 1) 



Uj—Ui—k 



A/ 



r(^, - + 1) 

+ qA^-'^^ (1 + 0(l/iV)) 



i + E 



CjkTjjyj - Ui + l) 
\ r(z/j - z/j - A; + 1) 



(25) 
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where Cj denotes a constant. The terms proportional to Cj arise from the regular part of 
the Borel transform and are not generally calculable. As expected r'f^'^ vanishes in the 
large N limit. 

Since mj^^ can be divergent under large limit, ~ m*^^^)j^-^rj^'* can be nonvan- 
ishing in large N limit, and thus we may write part of ( p2D as 

E[l-^^^Vrf) = 5:(Ai).,i^,+r,^ (26) 

3 j 

where we have separated those terms nonvanishing under large A^ limit from those vanish- 
ing, and put the former into and the latter into rj. Note that the nonvanishing 
part is linear in Ki. This is because the terms proportional to Cj in (^) give rise to a 
contribution only of 0{N~^^), thus affecting only r], which can be easily seen from the 
fact that (1 — m*^^))~-^ is at most as divergent as m!j^^ in the large A^ limit. Substituting 
(H) into (p2D we obtain 



We can now repeat this step by 



j j 
for a finite number of times (say / times) until I]j(l ~ '^Oij'^^j vanishes in the large A^ 
limit, to obtain 



K = lim 



(1 - Ai)-' ••■(!- Ai)-^(1 - m(^))-iK(^)l . (29) 



This is our main result. 

Now as an example, let us consider a Borel transform whose UV renormalon singularity 
is given by dH) with M = 2. Then from (|2|) and (|5]) 

m^f ^ = -— (1 + 0(1/A^)), otherwise mjf ^ = (30) 



and 



N""^-"^ (1 + 0{l/N)) + ciN-^^ (1 + 0{l/N)) 



T{V2 - Z/1 + 1 



rf) = -^^,N\l + 0{l/N)) + c,N-^-{l + 0{l/N)) (31) 
T[d + 1) 
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where 5 = z/i — z/2 — [z^i — z/2] — 1 < 0. Using the definition (^6[) we find 

/ n n \ 







(32) 



and obtain 









= hm 


(1 









;i-Ai)-vi 



m 



(Af)N-l 



/I 



(JV) 



(33) 



Now having given the series expression for the normahzation constants, a comment is 
in order. In deriving (|29|) we have assumed that all Ui > 0. As long as the normalization 
constants associated with positive Ui only are concerned, (p9D can be used without modi- 
fication. However, if some of z/j are negative, we have to consider A^^^b) instead of A(b) 
defined as 



M 



Y.K\(\-(5,b)-'^ 



i + $:c;,(i-/3o6)^' 



for 



\ ^34) 
Po 



such that "p satisfies t'j = z/j + p > for all z/j. Clearly Ki are related linearly with 
and so calculable once K[ are known. Now K[ from A^^^{k)) can be obtained in a series 
form using the steps taken in ([T7|) - (p9|), and is given by (|29|) with if — > if', Vi Ui and 
i(6) ^i(p)(6). 

In conclusion, we have shown that the normalization constants of the large order be- 
havior caused by an UV renormalon can be expressed, as in the case of an IR renormalon, 
in a calculable, convergent series, and thus can be computed to an arbitrary precision using 
the ordinary weak coupling expansion. Considering that the calculation of the normaliza- 
tion constants is equivalent to summing all sort of the higher order renormalon diagrams, 
it is surprising that they can be computed from the usual perturbation expansion. 
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